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Matrix elements of = 2 operators with Wilson fermions 
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We test the recent proposal of using the Ward identities to compute the Tf" — K'^ mixing amplitude with Wilson 
fermions, without the problem of spurious lattice subtractions. From our simulations, we observe no difference 
between results obtained with and without subtractions. From the standard study of the complete set of AS = 2 
operators, we quote the following (prehminary) resuhs: B^{2 GeV) = 0.70(10), (Or)^!,^^^ = 0.10(2)(l)GeV^ 
(08)<^=1 = 0.49(6)(0)GeV^ 



The main problem in lattice computations of 
the 4-fermion AF ~ 2 operators with Wilson 
fermions is related to spurious mixing among 
dimension-six operators. For example, chiral 
symmetry ensures that the operator responsible 
for the indirect CP- violation in the K'^ — sys- 
tem, Qi = s°7^(l - 75)^" s''7^(l - 75)^^, renor- 
malizes multiplicatively. Since the Wilson term 
explicitly breaks chiral symmetry, Qi mixes in- 
stead with all the other AS = 2 operators includ- 
ing 

Q,^-s''{l-l,)d<' s\l-i^)d'' 
Q,^r{l-^^)d'' s^(l-75)d" 

a =s'^(l-75)rf" s-''(l-f 75)^' 

=s"(l-75)rf'' s'(l + 75)d'' (1) 

where the superscripts denote color indices. The 
spurious mixing may seriously modify the chiral 
behavior of the operator Qi and hence need to be 
subtracted away, Q^^^^ \a) = Qi{a) + Ai{a)Qi{a) 
{i = 2, . . . , 5). After completing the subtraction 
procedure of the "wrong chirality" operators, the 
bare lattice regularized operator Q[^"'''\a) must 

*Talk given by Damir Becirevic. 



be multiplicatively renormalized (like in the con- 
tinuum): (5i(/i) = Zii{iJ,a)Q''i^''''' (a). A com- 
plete programme of renormalization of the oper- 
ators Qi-5 requires the computation of 9 renor- 
malization and 16 subtraction constants. This 
can be done perturbatively (see ref. Q]), but since 
the 0{as) terms are uncomfortably large, a non- 
perturbative determination of these 25 constants 
is mandatory. A theoretically simple method to 
implement the renormalization in the so-called 
(Landau)Rl/MOM scheme has been summarized 
in ref. In practice this programme is, however, 
quite complicated. Even more so since the pro- 
cedure included also the necessity to subtract the 
Goldstone boson (single and double pole) contri- 
butions ||-||] . That may cast doubts on the relia- 
bility of the method, i.e. on the level of control of 
the systematic errors. To address that issue, one 
needs a method allowing a computation without 
the necessity to subtract mixing with operators 
of the wrong chirality, and compare the results to 
the ones obtained by using the standard method 
(with subtractions). Recently, two such proposals 
appeared: twisted mass QCD |^ (see Q for the 
first numerical results), and method of the Ward 
identities B which we use is what follows. 
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1. Ward Identity Method ^ 

The method is extremely simple and to summa- 
rize it in a few lines we write the parity even (PE) 
operator Qi as Qi — VV + AA, in an obvious 
notation. For symmetry reasons, unlike the PE, 
the parity odd (PO) operators do not suffer from 
spurious mixings {i.e. Afj^ (a) = 0). It is then 
possible to apply the Ward identity on the matrix 
element of the PO operator, Qi = AV + VA, to 
get the matrix element of the PE one, Qi. By ap- 
plying the chiral rotation around the third axis in 
isospace, Su = 7511, Sd = —j^d {m^ — ma = to), 
on {P{x,tx)Qi{0)P{y,ty)), one arrives at 



gi(0) P{y,ty)) 



m{ 



J2 n(z,t.) P(f,t,) Qi(0) P{y,ty)) 



x,y,z,tz 



{J2P{x,t.) Qi(0) P{y,ty)) + (rot. sources) (2) 



where the terms arising from the rotation of the 
source operators P — ^755, vanish due to the 
symmetry under the charge conjugation. In the 
above identity, 11 = dj^d— uj^u. Thus, to get an 
informations on the r.h.s., we compute the l.h.s. 
where the operator Qi renormalizes only multi- 
plicatively. Moreover, Zp{fj,), which renormalizes 
the density 11, cancels against the one appearing 
in the mass renormalization constant, ZmifJ-) = 
Za/Zp{ij,), so that only Za is required. We com- 
puted the l.h.s. of the above identity on the lat- 
tice at /3 = 6.0 (16^x52) and at /3 = 6.2 (243x64). 
Results of both simulations were presented at the 
conference. Due to the lack of space, here we 
present only the results obtained on the finer lat- 
tice (/3 = 6.2, 200 configurations), where we work 
with Kg e {0.1339,0.1344,0.1349,0.1352}. Com- 
plete results with accompanying details will be 
presented in ref. 0. 

To confront the two methods, we studied the 
following ratios: 



7^l(^.) 



7^2(^.) = 



z\ {J2p{x, t,)P(o)) {J2 P{y, tv)Pm ' 



where, when needed, we used the non- 
perturbatively determined renormalization and 
subtraction constants. In the above formulae the 
time ty has been fixed, tx is left free and has 
been summed over all lattice. The role of denomi- 
nators is to eliminate the usual exponential terms 
in the numerator. 



O K=0.1339 
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Figure 1. Method without subtractions {TZi{tx)) is 
confronted to the standard one - with subtractions 
{'Ti2{tx))- Subtraction constants are computed non- 
pert urbatively. 

First, we observe that the plateaus for the ra- 
tio TZi{tx) exist and they are good for all the 
quark masses used in our study (see [||). Sec- 
ond, by using the method without (TZi(tx)) and 
with {TZ2{tx)) subtractions, we get the results 
very consistent with each other. In fig. |l| we 
show the ratio TZi/TZ2 for the smallest and the 
largest of our quark masses. From the observa- 
tion that (on the plateau) for all our masses we 
get 0.89 < 7^l/7^2 < 0.98, we conclude that the 
agreement of the two methods is satisfac- 
tory. Note that results of two methods have dif- 
ferent 0(a) effects. We checked that the chiral 
behavior for the matrix element {Ql{^J■)), as ex- 
tracted by using either of the two methods, is 
good, i.e. that TZ2{tx) when to^ ^ (see 
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fig. g). In addition, we performed the analysis of 
the data (w/o subtractions) when a small momen- 
tum is given to the external sources ("kaons"). 
After following the usual extrapolation proce- 
dures we get 

Bf^ = 0.70(10). (3) 

This error can be substantially reduced if the 
computation is made at several lattice spacing 
so that the momentum injections to the exter- 
nal kaons (for which the signals are noisier) are 
not needed ||]. 

2. S-parameters of the SUSY operators 

As seen in the previous section, the two meth- 
ods give very consistent results which makes us 
more confident that systematic errors introduced 
by using the standard method (with subtraction 
and renormalization constants computed non- 
perturbatively) , are indeed under control. We 
used the standard method to compute the matrix 
elements for the SUSY operators listed in eq. (p. 
They are parameterized as 



(7?°|Q2-5(/i)|if°) =S2-5(/i) C2_5X 

x(if"|s75d(M)0)(0|s75d(A*)|A'") 



(4) 



where C2 = —5/3, C3 = 1/3, C4 = 2, C5 = 2/3, so 
that the corresponding _B-parameters are unity 
in the vacuum saturation approximation. Each 
Bj-parameter is computed by replacing Q\—fQi 
in TZ2{tx) and dividing by Zp{ii) instead of Z\. 
Conversion from the RI/MOM scheme to the 
MS(NDR) scheme of ref. is made in pertur- 
bation theory with NLO accuracy. By linearly 
interpolating to the physical kaon mass, we ob- 
tain the following values: 

Sp(2 GeV) =0.64(5) (2), 



B. 



MS/ 



3_(2 GeV) =0.97(8)(12), 



.4 '(2 GeV) =0.87(6)(3), 



2 GeV) = 0.58(5)(1) 



(5) 



which, together with Bi given in eq. (j^), gives a 
complete set of _B-parameters of AS" = 2 opera- 
tors needed for the analysis of the SUSY effects 
in the — mixing amplitude 



3. Very briefly on ((7r7r)/=2 |g7,8(A^)|i^°) 

From the results of the previous section, one 
may get the useful information on the A/ = 3/2 
amplitude of the K — > tttt decay. After extrapo- 
lating to the chiral limit, the relevant bag param- 
eters are: B^ = Bf^{2GeY) = 0.46(6) (2), and 
B^ = Bp(2GeV) = 0.82 (7) (3). To predict the 
matrix elements in physical units, we will use the 
recipe of ref. [|l3| , which avoids the multiplication 
by the quark condensate (squared) and replaces it 
by the multiplication by Wp/^. Using that strat- 
egy, in the MS (NDR) scheme of ref. Q , we ob- 
tain: 

(77+107/^(2 GeV)|A:+) = -0.0193(33)(13) GeV^ 
(77+108^^(2 GcV)|A:+) = -0.092(10)(0) GeV^. 

These results, after using the soft pion theorem 
{i7r7r)j=2\Or,s\K°) = -(^+|0?i'|if+>/^^/2 @, 
lead to 

(07(2 GeV))^^'i = 0.104(18) (7) GeV^ 
(08(2 GeV))^'-'j, = 0.49(6)(0) GeV^, 

which can be compared to the 
preliminary estimates of the SPQcdR collabora- 
tion, (07(2 GeV))^lfJ^ = 0.021(11) GeV^ and 

(08(2 GeV))^lfj^ = 0.53(6) GeV^ as obtained 
by directly computing K ^ tttt matrix elements 
on the lattice. For comparison with other groups 
and other approaches, see ref. [T^ . 
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Figure 2. Chiral behavior of the matrix element of the operator ~ 1/a) as obtained by using the methods 

without (left figure) and with subtractions (right figure). The values of the fit parameters TZ = a + bX are also 
given {X oc niq, see ref. ||lO| for the precise definition). 



